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and a scalar field non-minimally coupled to matter. In the conventional approach to the stability 



We study the general properties of attractors in a cosmological model with tachyonic potential 



analysis the qualitative properties of the equations and of the long-term behavior of the solutions 
q , are investigated where the scaling solutions are the late time attractors and independent of the 

initial conditions. In this article, in a new approach, we examine the stability analysis of the model 
by simultaneously solving the dynamical system and best fitting the stability parameters with 
■ the observational data. The advantage of this approach is that the model which was of a purely 



mathematical nature become physically motivated. The number of the new critical points and 
also their properties, given in terms of the best fitted parameters, may alter due to the best fitting 



procedure. In a further step in stability analysis, we best-fit both the stability parameters and initial 
conditions with the observational data. The results impose more constraints on the trajectories in 
the phase space and provides more information about the dynamics of the cosmological model. 

X 

^ ; PACS numbers: 04.50.Kd; 98.80.-k 

Keywords: Tachyon, stability, phase space, distance modulus, statefinder, drift velocity. 



*Elcctronic address: hosscinf@guilan.ac.ir 
^Electronic address: a.salehi@guilan.ac.ir 



1. INTRODUCTION 



Recently, the observations of high redshift type la supernovae and the surveys of clusters 
of galaxies reveal the universe accelerating expansion and that the density of matter is very 
much less than the critical density . Whereas the above observational data properly 

complete each other, the dynamical dark energy (DE) proposal as an interesting possibility 
may arise to explain the observational constraints j^]. 

The two fine tuning and cosmic coincidence problems are the most serious issues with 
regards to the DE models and of the most frequently used approach to moderate these prob- 
lems is the tracker field DE scenario by employing scalar field models which exhibit scaling 
solutions (4]. Tracker models are independent of initial conditions used for field evolution 
but require the tuning of the slope of the scalar field potential. During the scaling regime, 
the scalar field energy density is of the same order of magnitude as the background energy 
density. The scaling solutions as dynamical attractors can considerably resolve the two 
above mentioned problems. By investigating the nature of scaling solutions, one can deter- 
mine whether such behavior is stable or just a transient feature and explore the asymptotic 



behavior of the scalar field potential Q . In studying the scalar field models with exponential 
potential there exist scaling attractor solutions js]. 

In quintessence dark energy model, there are two scaling solutions. One, the fluid-scalar 
field scaling solution, which remains subdominant for most of the cosmic evolution. It is 
necessary that the scalar field mimics the background energy density (radiation/ matter) in 
order to respect the nucleosynthesis constraint and can also solve the fine-tuning problem 
of initial conditions. Two, the scalar field dominated scaling solution, which is a late time 
attractor and gives rise to the accelerated expansion. Since the fluid-scalar field scaling 
solution is non-accelerating, we need an additional mechanism exit from the scaling regime 
so as to enter the scalar field dominated scaling solution at late times. For the discussion 
on the exiting —.one Canute 3. 

In scalar-tensor theories js| _ 12|, interaction of the scalar field with matter 18]- 21] and 



the presence of the tachyon potential in the formalism 2211- 231 separately are used to inter 



pret the late time acceleration. For tachyon dark energy 




as wel 



models in which a scalar field non-minimally coupled to the matter 



2n 



as other scalar-tensor 
, the scaling solutions 



have been investigated separately. In this paper, we integrate both non-minimal coupling 
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and tachyon models in one to investigate the scaling solutions and late time acceleration 
of the universe. The scalar field in our model is coupled with the matter lagrangian and 
also by its presence in tachyonic potential can be regarded as tachyon field. In our model, a 
prefect fluid with p m = 'ypm represent the matter in the universe, and the tachyonic scalar 
field candidates for DE. 

As discussed, in stability analysis we study the qualitative properties of the equations 
and of the long-term behavior of the solutions where the scaling solutions are the late time 
attractors and independent of the initial conditions. Here, in a new approach, we investigate 
the stability analysis of our model by simultaneously solving the system of equations and 
best fitting the stability parameters with the observational data. The advantage of this 
approach is that the model which was of a purely mathematical nature become physically 
motivated. The properties of the critical points that are now given exactly in terms of 
the best fitted parameters may change due to the best fitting effect. In a further step in 
stability analysis, we best-fit both the stability parameters and initial conditions with the 
observational data. The results impose more constraints on the trajectories in the phase 
space and provides more information about the dynamics of the universe. 

The well-known geometric variables, i.e. Hubble and deceleration parameters at the 
present time are used to explain the acceleration expansion of the universe. However, con- 
sidering the increased accuracy of the observational data during the last few years and 
generality of the DE models, new geometrical variables are introduced to differentiate these 
models and better fit the observational data. In this regard, a cosmological diagnostic pair 
{s,r}, called statefinder is introduced to differentiate the expansion dynamics with higher 
derivatives of the scale factor and is a natural next step beyond the well known geomet- 



ric variables 



29|. The statefinder pair has been used to explore a series of dark energy 



and cosmological models, including A cold dark matter (LCDM), quintessence, coupled 
quintessence, Chaplygin gas, holographic dark energy models, braneworld models, and so 



on 
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The Cosmological Redshift Drift (CRD) test which maps the expansion of the universe 
directly is also examined in here. We assume that the universe is homogeneous and isotropic 
at the cosmological scales 33| . The test is based on very simple and straightforward physics. 
Observationally, it is a very challenging task and requires technological breakthroughs, for 
more details see |34j-|38|. 
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The manuscript is organized as followers. In section 2, we derive the field equations for 
the model. In section 3, we perform stability analysis and obtain the autonomous equations 
in term of the new dynamical variables for the model. The critical points are also obtained, 
section 3.1 is designed to solve once again the autonomous equations by best fitting the 
stability parameters and initial conditions. Then we obtain the new properties for the 
critical points. The Stability of the best fitted critical points and phase space is studied 
in section 3.2. In section 4.1, for the best fitted model, the cosmological parameters are 
investigated. In addition, in section 4.2, the model is tested against observational data for 
the drift velocity. In section 5 we present summary and remarks. 



2. THE MODEL 



The model is presented by the action, 

R 



s = J - v(</>)y/i - + f(<f>)£ m ]V=gdx\ (l) 

where R is Ricci scalar, G is the newtonian constant gravity, and the second term in the 
action is tachyon potential. The f(<fi) is an analytic function of the scalar field. The last 
term in the lagrangian brings about the nonminimal interaction between the matter and the 
scalar field. The variation of action ([T]) with respect to the metric tensor components in a 
spatially flat FRW cosmology yields the field equations: 

W-^+JW (2) 
1-6 2 



2H + 3H 2 = - 1Pm f + V{<j>) A/ 1 - <j) 2 , (3) 



where we put 87rG = c = h = 1 and H — - with a is the scale factor of the universe, we 
also assume a perfect fluid with p m = 7p m . Note that in here 7 is the EoS parameter for 
the matter field in the universe. Variation of the action ([T]) with respect to the scalar field 
4> provides the wave equation for the scalar field as 



+ (1 - 2 )(3#0 +y) = - € -L{\ ~ 2 )l Pm , (4) 

where prime indicates differentiation with respect to and e = 1 — 37. From equations 
(EJ) (E]) on e arrives at the conservation equation, 

(p m f) + 3H(l + 1 ) Pm f = ep m f. (5) 



From equations fl2]) and by defining the effective energy density and pressure, p e // and 
Peff, one can identify an effective EoS parameter as 

Peff _ 7Pmf ~ V{4>)\Jl- (j? 

Peff Pmf 



— yeff _ irmj • VY-y V - r 

= 7— - -„ ; g7g lOj 



In the next section we study the stability analysis of the model in the phase space. 

3. PERTURBATION AND STABILITY ANALYSIS 

The structure of the dynamical system can be studied via phase plane analysis, by intro- 
ducing the following dimensionless variables, 

and parameters 

a=—. 8 = —. (8) 

Then using (|2])-(|5|), the equations for the new dynamical variables are, 

3 3t 3v 2 

x 1 = -2x(-- ~y x + JZ^) + eax ~ 3 ( X + ^ ( 9 ) 

y> = -2y(-- -^x+ + p y (fO) 

v yv 2 2 2 - 2x v ; 

Z > = -(1 - 2 2 )(32 + _ e (l _ ^2)|^ (11) 

Z Z 

where prime " ' "in here and from now on means derivative with respect to N = ln(a). By 
using the Fridmann constraint equation (j2J) in terms of the new dynamical variables, i.e. 

, + ^ = 1, (12) 



the equations (|9i)-(|TTl) reduce to, 



3 37 3y 2 

x = -2x(-- - —x + - _ — ) + eax - 3(1 + 7)2;, (13) 

v = -^-l-^+Y^ + ^y- ( 14 ) 



In term of the new dynamical variable we also have, 



* = _§ (1+7I+ ^,. (15) 



In stability formalism, by simultaneously solving x' = 0, y' = the fixed points (critical 
points) can be obtained. The critical points that depend on the cosmological and stability 
parameters 7, a and (3 are illustrated in Table I. 



TABLE I: The fixed points 



points FP1 


FP2 


FP3 


x 








y 







Substituting linear perturbations x' —> x'+5x', y' — > y'+8y', about the critical points into 
the two independent equations (U3l)- (jT4l) . to the first orders in the perturbations, gives us 
two eigenvalues Xi(i = 1, 2) which has to be negative as a requirement by stability method. 
In the following, the nature of the three critical points are given with the stability conditions: 



FP1 : X 1P1 = 3 + 0, X 2Pl = -3 7 + (1 - 3 7 )a 
FP2, 3 : A 1P2i3 = -6 - 2/3, A 2P2 , 3 = -3 - 3 7 + (1 - 3 7 )« - f3 



where 
FP1 : stablefor < 



a < 



a > 



FP2, 3 : stablefor < 



37 

1-37 
37 

1-37 
< -3 

37 



-3 
-3 



7<f 
1 



7>t 
1 = \ 



[[^ <a,/3> -3a 7 - 3 T + a - 3], [a < A, /3 > -3]] 7 < § 

:-3],[a> 3^T,/3 > -3]] 7 >i 



LL3 7 - 

/3> 



Y > a, (3 > — 3«7 — 37 
-3 7 = 1 



As can be seen, all the critical points are stable for the given conditions on a, (3 and 7. 
Ont the other hand, two of the cosmological parameters which relates the dynamics of 
the universe with the observational data are the EoS and deceleration parameters. In terms 
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of the new dynamical variables in our model they are given by, 

y 2 

U eff = JX+ ^—j, (16) 

q = ^ 1 + ^ + ^ri)- ( 1T ) 

Moreover, the statefinder parameters {r, s} in terms of the dynamical variables are, 

H d , H . , H . H , . 



r — 1 



(19) 



%-|) 

where in terms of the new dynamical variables is given by f[To"|) . For our model these 
parameters are presented in Table II: 



TABLE II: Properties of the fixed points 



points q u eff r s acceleration 



FP1 1/2 11 No 

FP2 -l-f-l-lf + M + i^ /3> _ 2 

FP3 — l — f — l — f ir + ir + i ~r /?>-2 



In the following, we solve the above equations for 7 = and 7 = 1/3 by best fitting the 
model parameters and initial conditions with the observational data using the x 2 method. 
This enables us to find the solutions for the above equations and conditions for the stability 
of the critical points that are physically more meaningful and observationally more favored. 
In the next section we best fit the model with the observational data for distance modulus. 



3.1. Best fitting the stability parameters and initial conditions 

The difference between the absolute and apparent luminosity of a distance object is given 
by, fi(z) = 25 + 5 log 10 d^z) where the Luminosity distance quantity, d L (z) is given by 

w-v+'^m- (20) 
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In our model, from numerical computation one can obtain H(z) which can be used to 
evaluate fi(z). To best fit the model for the parameters a and (3 and the initial conditions 
x(0), y(0), H(0) with the most recent observational data, the Type la supernovea (SNe la), 
we employe the x 2 method. We constrain the parameters including the initial conditions by 
minimizing the \ 2 function given as 

xL,.C, P, h {0)M o)) = E m 'f )M0)) ~ (21) 

1=1 1 

where the sum is over the SNe la sample. In relation (1211) . /i*' 16 and fi° bs are the distance 
modulus parameters obtained from our model and from observation, respectively, and a 
is the estimated error of the From numerical computation, Table III shows the best 

fitted model parameters for 7 = and 7 = 1/3. 

TABLE III: Best-fitted model parameters and initial conditions. 



parameters 


a j3 


x(0) y(0) £r(0) 




7 = 


-0.17 -1.62 


-0.2 -0.8 0.904 


548.8365694 


7 = 1/3 


-1.1 


-0.2 -0.8 0.909 


552.6788148 



Fig. 1) shows the constraints on the parameters a, (3 and the initial condition H(0) at 
the 68.3%, 95.4% and 99.7% confidence levels in both cases of 7 = 0, 1/3. Note that in case 
of 7 = 1/3, since equation ( JT4|) is independent of a, there is no constraint on a in this case, 
and therefore the equation becomes independent of this parameter as can be seen in the 
table III and Fig. 1. 
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a 





Fig. 1: The graph of confidence level for (a, H ) and H ) 
Left and Middle)for case 7 = 0, Right) for 7 = ~ 

Alternatively, we can plot the likelihood for the model parameters in both cases 7 
and 7 = 1/3 ( Figs. 2-4). 
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Fig. 2: The graph of 2-dim likelihood distribution for (a, H ) and (/3, H ) 
Left and Middle) for case 7 = 0, Right) for 7 = | 
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Fig. 3: The graph of 1-dim likelihood distribution for a, (3 and Hq for 7 = 



Fig. 4: The graph of 1-dim likelihood distribution for f3 and H for 7 = | 

The distance modulus, fi(z), plotted in Fig. 5, in both cases 7 = 0, 1/3 are best fitted 
with the SNe la observational data for the model parameters and initial conditions using \ 2 
method. 



V--z n~z 




0.5 1 1.5 0.5 1 1.5 



z z 

Fig. 5: The best-fitted distance modulus fi(z) plotted as function of redshift for 

left) 7 = right) 7 = | 

In the following we investigate the stability of the model with respect to the best fitted 
model parameters and initial conditions for the two specific choices of the EoS parameter 
for the matter in the universe, i.e. 7 = and 7 = 1/3. 

3.2. Stability of the best fitted critical points and phase space 

Solving the stability equations for the best fitted model parameters we find fixed points 
and the stability properties as illustrated in tables IV and V, for 7 = and 7 = 1/3 
respectively. 
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TABLE IV: Best-fitted fixed points for 7 = 



points (x, y) Stability 

FP1 (0, 0) unstable 

FP2 (0, 0.7) stable 

FP3 (0,-0.7) stable 

TABLE V: Best-fitted fixed points for 7 = \ 

points (x, y) Stability 

FP1 (0, 0) unstable 

FP2 (0, 0.8) stable 

FP3 (0,-0.8) stable 



From the above tables we see that, by best fitting the stability parameters, the critical 
points FP2 and FP3 are stable whereas FP1 is unstable in both 7 = 1/3,0. Thus, the 
number of critical points reduced from three to two. In addition, when best-fitting both the 
stability parameters and initial conditions, the best fitted trajectory for the system can be 
obtained which starts from an unstable critical point or infinity and approaches the stable 
critical point. In both graphs in Fig. 6), for 7 = 0, 1/3 all the best fitted trajectories for the 
stability parameters with different initial conditions leaving the unstable critical point PI in 
the past in the phase plane is shown moving towards the stable critical point P2 and P3 in 
the future. Moreover, the best fitted trajectories ( in red color) for both stability parameters 
and initial conditions are shown approaching the critical point P3. 
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Fig. 6: The phase plane for left) 7 = 0, right) 7 = | for best-fitted model parameter a and f3 

4. COSMOLOGICAL TEST 

To verify our model with the observation, in the following we perform two different kind 
of cosmological tests. 

4.1. Cosmological parameters 

In order to understand the behavior of the universe and its dynamics we need to study 
the cosmological parameters. We have best fitted our model with the current observational 
data by the distance modulus test. The cosmological parameters analytically and/or 
numerically have been investigated by many authors for variety of cosmological models. 
Applying stability analysis and simultaneously best fitting the model with the observational 
data gives us a better understanding of the critical points and the dynamics of these 
parameters. The effective EoS parameter, deceleration parameter and statefinders are 
given in terms of new dynamical variables in equations (I16til9j) . In tables VI and VII, the 
properties of the best-fitted critical points for 7 = 0, 1/3 are shown. 

In Fig. 7 and 8, the effective EoS parameters and deceleration parameter are shown in 
both cases, 7 = 0,1/3. All the trajectories shown in these graphs are best fitted for the 
stability parameters a and (3. In addition, the trajectory with red color are best fitted 
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for both stability parameters and also initial conditions. From Fig. 7 one observes that for 
7 = 0, 1/3 the phantom crossing do not occur in the past and future. In Fig. 7) left and table 
VI, for 7 = it shows that if both stability parameters and initial conditions are best-fitted, 
the universe starts from the unstable state at in the past with ou e ff = and approaches 
the stable state in future with u e ff = —0.47 (FP3). The best-fitted current value for the 
effective EoS parameter is about —5.5. 

From Fig. 7)right and table VII, for 7 = 1/3 we see that the universe start from unstable 
state in the past with u e ff — >■ —00 and approaches the stable state in future with u e ff = 
—0.43 (FP3). The best-fitted current value for the effective EoS parameter is about —6.2. 

From Fig. 8 for 7 = 0, we see that, the universe start from unstable state with positive 
deceleration parameter and approaches the accelerating stable state in future. It is 
interesting to see from the graph that for the best fitted stability parameters the universe is 
currently decelerating whereas for the best fitted stability parameters and initial conditions 
the universe is currently in accelerating state which is in compatible with the observational 
data. One see that from Fig.8 for 7 = 1/3, the universe starts from unstable state in the 
past with q — >• — 00 which is not compatible with observation. 



TABLE VI: Properties of the best-fitted fixed points for 7 = 
points q oj e ff r s acceleration 

FP1 1/2 11 No 
FP2 -0.2 -0.47 -0.11 0.54 yes 
FP3 -0.2 -0.47 -0.11 0.54 yes 
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Fig. 7: The graph of best-fitted effective EoS parameter, oo e ff for left) 7 = 0, right) 7 
The red trajectory is for the both best-fitted model parameters and I.C.s 
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Fig. 8: The graph of best-fitted deceleration parameter, q for left) 7 = 0, right) 7 = |. 
The red trajectory is for the both best-fitted model parameters and best I.C.s 

Figs. 9-13 show the statefinder diagrams {s, q} and {r, q} and their evolutionary trajecto- 
ries. As above, in all these graphs, the statefinder trajectories are best-fitted for the stability 
parameters, whereas the trajectories with red color are also best fitted for initial conditions. 

From the graph of the statefinder {r, q}, {s,q} and {r,s} for 7 = we see that all the 



TABLE VII: Properties of the best-fitted fixed points for 7=3 
points q w e// r s acceleration 

FP1 1/2 11 No 
FP2 -0.45 -0.63 -0.02 0.37 yes 
FP3 -0.45 -0.63 -0.02 0.37 yes 
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best-fitted trajectories start from the standard cold dark energy (SCDM) in the past which 
is an unstable critical point and tend their evolution to the stable state in the future. The 
current values of the trajectories are also shown and can be compared with the position of 
SCDM. However the result is not promising for the case 7 = 1/3 as the trajectories start 
from not a physically known state in the past. 




-0.4 -0.2 0.2 0.4 0.6 -0.5 -0.4 -0.3 -0.2 -0.1 



q q 
Fig. 9: The graph of the best-fitted statefinder {s, q} for left) 7 = 0, right) 7=3- 



The red trajectory is for the both best-fitted model parameters and best I.C.s 




Fig. 10: The graph of the best-fitted statefinder {r, q} for left) 7 = 0, right) 7 = |. 
The red trajectory is for the both best-fitted model parameters and best I.C.s 
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Fig. 11: The graph of the best-fitted statefinder {r, s} for left) 7 = 0, right) 7 = - 
The red trajectory is for the both best-fitted model parameters and best I.C.s 
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Fig. 12: The graph of the best-fitted statefinder (r, ln(a)) for left) 7 = 0, right) 7 = |. 
The red trajectory is for the both best-fitted model parameters and best I.C.s 
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Fig. 13: The graph of the best-fitted statefinder (s,ln(a)) for left) 7 = 0, right) 7 = |. 
The red trajectory is for the both best-fitted model parameters and best I.C.s 
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4.2. CRD test 



The Cosmological Redshift Drift is extracted from 



1 + z)H - H(z) 



(22) 



known as Mc Vittie equation. The equation immediately leads to velocity drift 

cH{z) 



v = cHq 



l + z 



(23) 



with respect to the redshift that, by using H(z) from numerical computation in our model, 
can be obtained against observational data. In our model, H(z) is taken from numerical 
calculation in terms of the stability dynamical variables and after best fitting the model 
with the stability parameters and initial conditions. Fig. 14 shows a comparison between 
the velocity drift in our model for 7 = 0, 1/3 and ACDM model. As can be seen our model 
for 7 = better fits the observational data in comparison with ACDM mode. It also shows 
another disadvantage in 7 = 1/3 case. 




ACDM 

Best fit of SNIafory=0 
Best fit of SNIa for j=- 



Fig. 14: The graph of velocity drift as function of redshift for 7 = and 7 



5. SUMMARY AND REAMRKS 

This paper is designed to study the attractor solutions of a cosmological model with tachy- 
onic potential and non minimally coupled scalar field with the matter, by stability analysis 
and making use of the 3-dimensional phase space of the theory. The model characterized 
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by the scalar field <fi , the scalar potential V ((/)), and the scalar function f((p) nonminimally 
coupled to the matter lagrangian. The scalar function f(<j)) and potential V(4>) are arbitrary 
and not specified in advance. The stability analysis gives the corresponding conditions for 
tracking attractor and determines the type of the universe behavior in the past and future. 

In this work, with a new approach in stability analysis, we simultaneously best-fitting the 
stability parameters with the observational data using \ 2 method and solving the system 
of equations. The advantage of this approach is that the obtained critical points are obser- 
vationally verified and thus physically more promising. The best fitted critical points are 
physical critical points. By this approach in our model we find that one of the mathemat- 
ically known stable critical point (FP1) is now physically unstable. Another advantage of 
this approach is that the fixed points coordinates are best fitted. We added even one more 
step in the process and best-fitted both stability parameters and initial conditions in the 
model. As a result, we obtained the observationally verified trajectories in the phase plane. 

We then study the cosmological parameters such as effective EoS parameter, u e ff, decel- 
eration parameter, q, and statefinder parameters for the model in terms of the best-fitted 
stability parameters. It shows that the EoS parameter for both 7 = 0, 1/3, does not cross the 
cosmological divide line in the past and future. In 7 = case, the universe starts from un- 
stable state in the past with u e ff = and finally tends to the stable state in the future with 
u e ff = —0.47. The current value of the EoS parameter of the universe for both best fitted 
stability parameters and initial conditions is within the range of observationally accepted 
values. The model also tested with CRD observational data and shows a better match with 
the data compare to the ACDM model and also case 7 = 1/3. With both best-fitted stability 
parameters and initial conditions, the deceleration parameter q satisfies q < at present and 
q > in the matter dominated era, in the case of 7 = 0. The statefinder parameters show 
that the universe start from SCDM in the past and approaches a state near SS in future for 
7 = 0. In 7 = 1/3 case, the model is not observationally verified as the universe starts from 
an accelerated state in the past and the EoS parameter also begins from an unstable state 
with u e ff = —00. In addition, the model does not explain the CRD data. Note that in 
this case, from stability equations f lTBUT^ . the new dynamical variables are independent of 
one of the stability parameters, a. Since a is proportional to 4, the scenario is similar to a 
cosmological model with no coupling between the scalar field and matter which can justify 
the unsatisfactory result in 7 = 1/3 case. 
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